o 
o 

(N 

1—5 

(N 






Emptiness formation probability at finite temperature for the 

isotropic Heisenberg chain 

Frank Gohmann, * Andreas Kliimper and Alexander Seel 

Bergische Universitdt Wuppertal, ^2097 Wuppertal, Germany 



-)— > 

a 

o 
o 



Abstract 

We present an integral formula for a special correlation function of the isotropic spin- 1/2 antiferromagnetic 
Heisenberg chain. The correlation function describes the probability for the occurrence of a string of consecutive 
up-spins as a function of temperature, magnetic field and length of the string. 
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In a recent work [3] we derived (multiple) integral 
representations for a generating function of the a^-a^ 
correlation functions of the spin-i XXZ chain, 



_iaJ + <_,aJ + A(a|_ia|-l) , (1) 



at finite temperatures T and finite values of an exter- 
nal, longitudinal magnetic field h. Our formulae gen- 
eralize the zero-temperature results of [7] . 

In order to relate our work to experiments the in- 
tegrals should be calculated. Although this task may 
appear rather challenging at first sight, the recent 
progress made in the calculation of related integrals 
(see e.g. [1,2,13]) and their asymptotic behaviour [6,12] 
at zero temperature lets us hope that we eventually 
will be able to extract numbers from our formulae. 

The integrals [4,5,11] calculated in [1,2,13] do not 
describe local correlation functions but density matrix 
elements (termed 'elementary blocks' in the mathemat- 
ical literature) . One of these non-local correlation func- 
tions, called the emptiness formation probability [11], 
is contained as a special case in our formula obtained 
in [3]. In [3] (see eqs. (125), (126)) we represented 
it as a multiple integral with symmetric integrand of 



similar form as in [7]. Below we obtain another inte- 
gral representation with asymmetric integrand which, 
when further specialized to the isotropic antiferromag- 
net (J > and A = 1 in (1)), is the finite tempera- 
ture generalization of the formula obtained in [11] that 
served as a starting point for the explicit calculation of 
the emptiness formation probability in [1]. We there- 
fore believe that it may also be the appropriate starting 
point for proceeding further to finite temperatures. 



1. Thermodynamics 

There are several alternative possibilities of repre- 
senting the free energy of the Heisenberg chain (1) as 
an integral over auxiliary functions. The representation 
based on the quantum transfer matrix [14, 15] and de- 
rived in [9, 10] turns out to be most appropriate in the 
present context, since the auxiliary functions involved 
also appear quite naturally in the derivation of multi- 
ple integral representations of correlation functions [3] . 
Using the formulation of [3] and specializing to J > 
and A = 1 we can write the free energy per lattice site 
in the thermodynamic limit as 
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fih,T) 



T 



duj ln(l -|- o(a;)) 
27r uj(lu + i) 



(2) 
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Fig. 1. The canonical contour C 

The contour C encloses the real axis at a distance 
I7I < |. The auxiliary function o(A) is the solution of 
the non-linear integral equation 



h 



lna(A) = -- 



2J 



duj ln(l + a(<^)) 



T A(A + i)T y^ '^ 1 + (A -6^)2 



(3) 



2. Emptiness formation probability 

The emptiness formation probability [11] is defined 
as the thermal average 



nm) = (nr=i(JE! 



3'T- 



(4) 



Within the formalism developed in [3] it is obtained 
as the homogeneous limit 



P(m) 



ei,...,e 



lim P(m|te}™ 



j^iJ 



(5) 



of a more abstract correlation function P(m|{5j}™^i). 
Here we present a new integral representation for this 
function, 



P{m\{^,}T=i) 



dcui 



27r(l + a(a;,))Jni 



det G{lj J, ^k) 



m 

[n 



<j<,k<m 



(6-G) 



El 



<j<,k<m 



{ujj —UJk+l) 



(6) 



As in our general formula in [3] the function G(A, ^) 
appearing in the determinant on the right hand side 
has to be calculated from a linear integral equation 
which in the isotropic case (A = 1) reads 



G(A,0 = - 



(A-C)(A-C-i) 
duj 1 



+ 



G{lu,0 



c TT 1 + (A -6^)2 l + a{ui) 



(7) 



The formula (6) follows from equation (111) of [3] by 
(i) taking the limit ip — > — cx3, (ii) using appendix C 
of [7], and (iii) taking the isotropic limit. 



Taking the homogeneous limit of (6) is not diffi- 
cult [8]. We obtain the expression 



p{m) = [n 



duji 



l\Jc2n{l + a{iJ,)) 



■det 



4 'G(^„e)i,^o nr=i(^-i)'"'^. 



— 1 m — i 

'■' . 
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(fe-1)! ni<i<fe<™(wi -Wfe+i) 



(8) 



This is the direct generalization to finite temperatures 
of the formula obtained in [4, 8, 11]. 



3. Reduction of the integrals 

We believe that either (6) or (8) or the symmetrized 
form of (6) is a good starting point for studying the 
possibility to reduce the multiple integrals to simpler 
ones. That this may be possible in principle can be in- 
ferred from the example m = 2, where, for h = Q and 
in the isotropic limit, the emptiness formation proba- 
bility P(2) can be expressed in terms of the internal 

energy e, 

1 

(9) 



^(2)4+if7- 



Inserting (2) into the thermodynamic expression for 
the internal energy e — 9i/t(//T) we obtain 



P(2) 



1 
12 



did 



G(CJ,0) 



(o; — i) 1 -I- o((j;) 



(10) 



Note that there is only a single integral on the right 
hand side. Specializing (8) for m = 2, on the other 
hand, yields 



P(2) 



duj\ 



duJ2 



a;i(a;2 



47r2 Jq 1 + a{uji) Jc 1 + o-i^'i) uJi — u;2 + i 



(11) 



which involves two integrations. 



4. Discussion 

The integral representations (6), (8) are likely can- 
didates for extraxting for the first time explicit num- 
bers for finite temperature next-to-nearest neighbour 
correlation functions of the isotropic Heisenberg chain 
and may serve as a model case to initiate a systematic 
study of the more general integrals derived in [3] . 
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